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ABSTRACT 

We have developed a semi-analytic method of parameterizing iV-body simulations 
of self-gravity wakes in Saturn's rings, describing their photometric properties by means 
of only six numbers: three optical depths and three weighting factors. These numbers 
are obtained by fitting a sum of three gaussians to the results of a density-estimation 
procedure that finds the frequencies of various values of local density within a simulated 
ring patch. 

Application of our parameterization to a suite of iV-body simulations implies that 
rings dominated by self-gravity wakes appear to be mostly empty space, with more than 
half of their surface area taken up by local optical depths around 0.01. Such regions 
will be photometrically inactive for all viewing geometries. While this result might be 
affected by our use of identically-sized particles, we believe the general result that the 
distribution of local optical depths is trimodal, rather than bimodal as previous authors 
have assumed, is robust. 

The implications of this result for the analysis of occultation data are more concep- 
tual than practical, as we find that occultations can only distinguish between bimodal 
and trimodal models at very low opening angles. Thus, the only adjustment needed 
in existing analyses of occultation data is that the model parameter r gap should be in- 
terpreted as representing the area-weighted average optical depth within the gaps (or 
inter- wake regions), keeping in mind the possibility that the optical depth within those 
inter-wake regions may vary significantly. 

The most significant consequence of our results applies to the question of why "pro- 
peller" structures observed in the mid- A ring are seen as relative-bright features, even 
though the most prominent features of simulated propellers are regions of relatively 
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low density. Our parameterization of self-gravity wakes lends preliminary quantitative 
support to the hypothesis that propellers would be bright if they involve a local and 
temporary disruption of self-gravity wakes. Even though the overall local density is 
lower within the propeller-shaped structure surrounding an embedded central moonlet, 
disruption of the wakes would flood these same regions with more "photometrically ac- 
tive" material (i.e., material that can contribute to the rings' local optical depth), raising 
their apparent brightnesses in agreement with observations. We find for a wide range 
of input parameters that this mechanism indeed can plausibly make propellers brighter 
than the wake-dominated background, though it is also possible for propellers to blend 
in with the background or even to remain dark. We suggest that this mechanism be 
tested by future detailed numerical models. 

Subject headings: Saturn, rings 

Running header: Analytic parameterization of self-gravity wakes 



1. Introduction 



Saturn's dense A and B rings are pervaded by a microstructure, dubbed "self-gravity wakes," 
of alternating dense and rarified regions that arise due to a rough balance between the clumping 
together of particles under their mutual self-gravity and their shearing apart again due to tidal 



forces (Julian and Toomre 1966 Salo 1992; Richardson 1994 1 . The characteristic alignment of 



elongated self-gravity wake structures causes variation with viewing geometry in the brightness of 
the rings as seen in ground-based optical (Camichel 1958 Salo et aL"||2004 French et al.||2007 ) and 
microwave observations ( Dunn et al.||2002 2004 1 as well as spacecraft images ( Franklin et al.fl987 
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properties. 



Existing analyses of stellar occultation data from Cassini UVIS ( Colwell et al.|20"06 2007) and 
Cassini VIMS ( Hedman et al.|2007 Nicholson and Hedman|2009 1 explain the observations through 
the use of simple models that assume an optical-depth dichotomy, with nearly-opaque wakes (with 
optical depth T wa k e ) and a low but relatively constant optical depth in the spaces between the wakes 
( T gap)- However, while some preliminary work has shown that such a model can produce brightness 
profiles that are consistent with data ( Hedman et al.|2007 Salo et al.|2008 ), it has not been verified 
that such a bi-modal model describes the actual nature of simulated wakes, not to mention real 
ones. And if it does, how do r wa k e and T gap relate to environmental parameters such as overall 
occultation optical depth and coefficient of restitution? What do the values of T gap inferred from 
observations tell us about the conditions under which wakes occur? 



Salo et al. (2004) and Porco et al. (2008), on the other hand, produce plots of ring brightness 
at various geometries for various simulation input parameters, inferring the best input parameters 
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by finding simulated brightnesses that best compare to observations. But where is the brightness 
coming from? For images in reflected light, is the brightness varied primarily by changes in the 
opacity of dense wake structures, or by the fractional area covered by them? For images in trans- 
mitted light, do more photons get scattered into the camera from the rarefied regions between the 
wakes (the wakes themselves being largely too opaque to allow much light to pass through them), 
or are more photons being scattered by the edges of the wake structures? Better understanding of 
these questions would help to guide future simulations, and might allow one to infer more from the 
available observational data. 

Yet another current question concerns the nature of "propeller" structures observed in images 



of the mid-A ring (Tiscareno et al. 2006, 2008; Sremcevic et al. 2007), each of which is inferred to 



be a disturbance surrounding an unseen moonlet embedded within the ring. Observed propellers 
generally show the disturbed region as brighter than the background ring brightness, though theo- 
retical models ( Spahn and Sremcevic|200"0 Sremcevic et al.|2002 Seifi et al.|2005 Lewis and Stewart 



2009) indicate that the most prominent feature of the disturbance is a decrease in local density. 



Among other hypotheses for this curious behavior, Tiscareno et al. (2006 2008) have proposed 



that self-gravity wakes tend to lock up ring material into a photometrically inactive state, and that 
propellers can release this material by locally and temporarily disrupting the wake structure, thus 
causing the disturbed area to contain more photometrically active material even if it contains less 
material overall. Our parameterization of self-gravity wakes allows us to evaluate this hypothesis 
in a more quantitative fashion than has previously been possible. 

A major reason why these questions remain largely unanswered lies in the difficulty of local 
density estimation, the quantification of an underlying density function based on the locations 
of a finite number of particles, for a strongly heterogeneous density distribution. The simplest 
method of density estimation is to divide a simulation space into bins of equal size and to count 
the number of particles in each bin, but this fails to give meaningful results for current simulations 
of self-gravity wakes. This is because the rarefied regions require large bins in order to accumulate 
enough particles in each bin to give good statistics, but such large bins will badly smear the sharp 
boundaries between the dense wakes and the rarefied regions. Contrariwise, if the bins are small 
enough to resolve the sharp boundaries, nearly all bins in the rarefied regions will contain zero 
particles (with a few containing one particle), giving no sense of the average density in the rarefied 
regions. 

We have constructed a density estimation method, to be applied to self-gravity wake simula- 
tions, based on circular bins that expand in rarefied regions and contract in dense regions. A given 
circular bin is only used if the particles within it satisfy certain criteria (described in detail below) 
to ensure that they are sufficiently evenly distributed; if those criteria are met, then the bin grows 
until they are no longer met. The overlapping circular bins are then projected onto a finely-meshed 
rectangular grid; nearly all bins in the latter grid contain multiple overlapping circular bins, and 
the value of the density at each location is taken from the average of the overlapping circular bins. 
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Applying our density estimation method to a set of numerical simulations of a ring patch 
characterized by self-gravity wakes, we parameterize the general properties of each simulation as 
a weighted combination of three optical depth values. This semi- analytic treatment allows us to 
comment on trends that are applicable to a wide range of input conditions. We specifically discuss 
the implications of our results for previous interpretations of occultation data, as well as for the 
question of why "propellers" appear as relative-bright features. 

Section [2] describes our suite of numerical simulations, the results of which constitute the input 
data for our semi-analytic method. Section [3] describes our density estimation method. Section [4] 
describes our results, and Section [5] provides further discussion. Section [6] presents a summary and 
conclusions. 



2. Wake Simulations 

We carried out a series of numerical simulations of an orbiting patch of ring particles. Details 



of the numerical technique are provided in Porco et al. (2008), but are briefly summarized here 



for convenience. Particle trajectories are computed in a rotating frame (the patch) using Hill's 
equations of motion with self-gravity. The orientation of the patch is such that the x-direction 
points radially away from the planet, the y-direction is in the direction of the patch center motion, 
and z is perpendicular to x and y according to the right-hand rule. Boundary conditions are applied 
in the x and y directions to keep the particle number constant. Duplicates of the patch surround 
it in the orbital plane to provide a smoother gravity potential and to allow for collisions at the 
patch boundary. Duplicates in the ±x directions are offset in ±y to account for differential shear 



(Wisdom and Tremaine 1988). A second-order leapfrog integrator adapted to the rotating frame 
is used to integrate the equations of motion. Collisions are predicted during the "drift" update 
of each integration step and are carried out using billiard-ball restitution equations. We ignored 
surface friction and particle spin. We chose a timestep of 5 seconds, which is more than 1,000 times 
smaller than either the orbital period or the dynamical interaction time between two particles, 
wl/Gp. Initial conditions consist of a uniform distribution of dynamically cold particles in a thin 



Table 1: Parameters for our simulations of self-gravity wakes. 
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Tdyn is input mean dynamical optical depth, a is input mean surface density, 

N is total number of particles, and A cr = 47t 2 G(t/k 2 is the Toomre critical wavelength 
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slab; equilibrium, measured as a flattening of the components of the velocity dispersion in the 
patch, is established typically within 10 orbits of the central patch around Saturn. 

We used seven values of the input mean dynamical optical depth r^ yn , defined as the total cross- 
section area of particles divided by the total patch area (Table [TJ. All simulations used a Saturn- 
centered orbital distance of 130,000 km, a monodisperse particle-size distribution (i.e., identically- 
sized particles) of radius R = 1.667 m, and internal density equal to the local Roche critical density 
0.45 g cm -3 ( Porco et al.|2007 20081. For this combination of particle properties, which are typical 
values for the mid-A ring, the numerical values for the input mean surface densities (in units of 
g cm -2 ) are simply 100 times the numerical values for r^ yn . The patch dimensions of 510 x 1260 m, 
with the short axis oriented in the radial direction and the long axis along the orbital direction, were 
chosen to be always greater than 4 x 10 times the Toomre critical wavelength (A cr = Ait 2 Gg/ re 2 , 
where G is Newton's constant and n is the epicyclic frequency), even for the simulations with 
the highest densities, which ensures that no individual particle or structure can reach across the 
periodic boundary conditions to meaningfully interact with itself. 



Other investigators have used a parameter r p (Ohtsuki 1993 Salo 1995) or r* h (Daisaka et al 



2001), a modified ratio of Hill radius and particle radius, to track the susceptibility to wake for- 



mation of a given simulated ring patch. This parameter depends only on the particle's internal 
density and semimajor axis (Daisaka et al. 2001 Eq. 8). Furthermore, the Roche critical density, 



at which a particle fills its own Hill sphere, depends only on the particle's semimajor axis (Porco 



et al.||2007 Eqs. 1 and 2). Combining these two, any particle at the Roche critical density will have 
; ~ 0.85. One can now locate our suite of simulations as a vertical line in the two-dimensional 



parameter space shown in Fig. 2 of Daisaka et al. (2001 ), and can calculate from their Eq. 13 that 



the onset of self-gravity wakes should occur for r > 0.25. 



For the primary batch of simulations, we used the coefficient of restitution stated by Borderies 
et al. (1984), building on the work of Andrews (1930), with v* = 0.001 cm s . This relatively 



dissipative law, which yields lower post-collision speeds than the more commonly-used coefficient- 



of-restitution law formulated by Bridges et al. (1984), was identified in the analysis of Porco et al. 



(2008) as leading to a better fit to observations of this particular region of the A ring. Furthermore, 



we carried out one additional simulation that was identical to the one described above with Td yn = 
0.5 and a = 50 g cm~ 2 , but using the Bridges et al. (1984) law. The latter simulation is intended 



to be identical to that used by Salo et al. (2004) for their photometric analysis. Throughout this 



paper, comparison of the results given by the two coefficients of restitution can be thought of as 
giving a crude approximation of the variation that might arise in our models due to uncertainty of 
various input parameters. 



1 2 3 4 5 
Particle # 



2 4 6 8 10 
Particle # 



Fig. 1. — Illustration of our density estimation method using a sample snapshot from the simulation 
with a = 50 g cm~ 2 at t = 30 orbits (top left), at a location centered at xo = 84 meters and 
yo = 47 meters from the frame center (inset at top right). The first circle (green circle, including 
particles shown in green) was validated by a "center of mass" less than 0.5 m from the micro-bin 
center (i.e., the green cross is inside the blue circle) and also by (lower left-hand plots) radial 
and azimuthal distribution of particles with standard deviations not too far from uniformity. The 
attempt to expand the circle (red circle, including particles shown in red and green) failed due to 
a "center of mass" too far from the micro-bin center (i.e., the red cross is outside the blue circle), 
and also due to (lower right-hand plots) excessive deviation from uniformity in both the radial and 
azimuthal distribution of particles. In the plots, circles indicate the distribution of r 2 values and 
squares indicate values. 
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3. Density Estimation Method 

After each simulation ran for a long-enough period to reach equilibrium (i.e., all components 
of the velocity dispersion have settled to roughly constant values), we took seven snapshots at 
t = {18, 20, 22, 24, 26, 28, 30} orbits. For each snapshot, containing the positions of all particles at a 
given time, we performed the following algorithm on a central sub-patcrj^jof dimensions 300 x 750 m. 

We first projected the three-dimensional particle positions onto a plane. We used the ringplane 
{z = 0) for this purpose, though this need not be the case. In order to quantify the distribution of 
densities observed at a slant-path through the ring, any plane perpendicular to a desired line-of-sight 
may be used. 

We divided the sub-patch into "micro-bins" of size 0.25 m, small enough to clearly resolve the 
boundaries between wake and inter- wake regions, forming a grid with dimensions 1200 x 3000 micro- 
bins. At each micro-bin location [xo,yo]j we begin by drawing a circle around the five particles 
closest to that location. The radius of the circle we denote as rc, and each of the N particles within 
the circle is assigned a position [r, 0] relative to the center of the circle and the x = axis. We then 
perform three tests to determine whether the particles are evenly distributed (Fig. [T]) : 

1 . The average position of the particles inside the circle (the "center of mass"[^]) must be within 
0.5 meters of the center of the micro-bin (the blue circle in Fig. [T|; 

2. The number of particles whose distance from the center of the circle is less than or equal to r 
should increase linearly with r 2 . That is, for a perfectly even distribution, the expected value 
of the radial location of the nth particle would be r n = r c {n — 1/2) /N, where n ranges from 1 
to N. We compute the standard deviation of the residual between the actual and ideal values 
of r n , which must be less than the total range divided by the number of particles (tq/N) to 
pass the test. 

3. Similarly, for a perfectly even distribution, the angular separation between adjacent particles 
would be constant, which is to say that the expected value of the longitude of the nth particle 
would be n = 2ir(n — 1/2) /N. We find the standard deviation of the residual between the 
actual and ideal values of n , which must be less than the total range divided by the number 
of particles (2ir/N) to pass the test. 

If all three tests are passed, then the micro-bin is validated and used in the final solution. 
Furthermore, we try to expand the circle, increasing rc to take in the next five particles and 
redoing the above-mentioned tests on the new distribution. If the tests are passed again, the 



*Our reason for focusing on a sub-patch was simply to reduce computation time. 

'This is a figure of speech; in the case of a particle-size distribution, when the average position of particles and the actual 
center of mass may not be the same, it is particle locations with which we are concerned. 
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iteration continues until the maximum value of tq = 12.5 meters is reached. If the tests are not 
passed, then we revert back to the last successful value of rc for that location. 

We identify five parameters in our method that require human input. These include the 
number of particles added to the circle in each iteration (we chose five, which is large enough for 
each iteration to have reasonable statistics, but small enough to allow the validation of bins that 
are close to a wake/inter- wake boundary), the maximum value of rc (we chose 12.5 meters, which is 
less than the characteristic spacing between the dense wake structures), and the threshold criteria 
for the three tests enumerated above (the values described could reasonably be multiplied by a 
scalar of order unity). Aside from this human input, our method is fully automated. 

This process is illustrated in Fig. [T] where the first circle (green) results in a successful valida- 
tion, while the attempt to expand the circle (red) fails. In the top-right panel, Test #1 can be seen 
in that the average position of the first five particles (the green cross) is inside the blue circle, while 
the average position of the first ten particles (the red cross) is not. Tests #2 and ^3 are illustrated 
in the plots, where the horizontal dotted lines indicate the standard deviation criteria. Note that 
the plotted residuals for r 2 and 9 have been normalized by the threshold criterion values; when the 
circle is expanded, the threshold criterion for 6 decreases simply because TV" is twice as large, while 
the threshold criterion for r 2 slightly increases because a larger r c compensates for the increased 
N. 

In the gaps between the dense wake structures, it commonly happens that even the maximum- 
size circle (rc = 12.5 meters) contains fewer than five particles. In this case, we designate the 
location as "sparse" and proceed with the particles that do occur in the maximum-size circle. For 
"sparse" bins, we waive Test #1 (the "center-of-mass" test), in order to avoid large areas of the 
inter- wake gaps having no coverage at all with valid bins. 

When the above process is complete, each location on our 1200 x 3000 grid has one of two 
states: a valid circle of radius rc with a certain number of particles within it (which can be resolved 
into a density of number of particles per unit area), or no valid result. We then visualize the valid 
circles and find that nearly every location within the grid is covered by one or more overlapping 
circles. Finally, we calculate the density at each location on the grid by summing the particles in 
all circles overlapping that location and dividing by the sum of the areas of the same overlapping 
circles 



Fig. 2. — Sample result of our density estimation method, showing a portion of the same snapshot 
shown in Fig. [I] Red dots are particle locations, which constitute the inputs for density estimation. 
Calculated densities are shown as a grayscale background ranging from black (low density) to white 
(high density). 



4. Results 

A sample result of our density estimation method is shown in Fig. |2j We find that regions 
with high calculated density (bright in the figure) correlate well with regions with tightly-clustered 
particles, and that regions with low calculated density (dark in the figure) correlate well with 
regions in which particles are sparse. Although some instances of intermediate calculated density 
(gray in the figure) are on the boundaries between dense and sparse regions, a smearing effect that 
our method was designed to avoid, we find that the large majority of regions with intermediate 
calculated density genuinely have particles clustered to an intermediate degree. 

Now that we have a grid of local densities within a self- gravitating ring patch, we can make a 



' Rather than simply averaging the densities in the overlapping circles, this method weights the result somewhat towards 
those circles that have more particles, and thus better statistics. 
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Fig. 3. — Histograms of local dynamical optical depth, Td yn (x, y), for the eight simulations described 
in Section [2j calculated using the method described in Section [3] Least-squares fits are made to 
a sum of three gaussians (or, in the top-left case, two gaussians). Within each plot are listed the 
fit parameters: the height, center location, and full-width at half-maximum (the latter in units of 
log r) of each gaussian. 



histogram showing the frequency with which each value occurs (Fig. [3]). For this purpose, we express 
the local density in terms of the local dynamical optical depth, r^ n (x,y) = p(x,y)7rR 2 /£ 2 , where 
p(x,y) is the number density (particles per bin) calculated by the method described in Section |3j 
i = 0.25 m is the bin size, and R = 1.667 m is the particle radius. Note that this is not the same 
quantity as that measured in optical observations, the photometric optical depth T p hot = 

-InT, 

where the transparency T is the fractional area not blocked by particles. The two are nearly equal 
at low values, but diverge at large optical depths when particles come together closely enough that 
their inability to occupy the same location in space becomes important. 

Some of the histograms in Fig.|3]have an anomalously high value in the left-most bin, containing 
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the lowest optical depths. This is due to a small fraction of bins which, in the application of our 
density estimation method, were not overlapped by any valid circular bins. This generally occurs 
in small regions of very low density that are surrounded by regions of higher density; such regions 
are visible in Fig. [2] as uniformly black. Rather than further refine the method to better account 
for these difficult cases, we simply assigned to these locations a density equal to the lowest densities 
obtained. The fact that these left-most spikes are not very large in any of the histograms justifies 
our decision to neglect this effect as a small perturbation. 

We obtain quite good results by modeling the histograms in Fig. [3] (using log Td yn as the 
independent variable) as a sum of two or three gaussians. We justify the choice of a gaussian fit as 
follows: The results of our density estimation, though quite good, do contain some variation within 
regions that should have a single density; sometimes the local density is higher in the immediate 
vicinity of a single particle, and sometimes the density decreases smoothly with distance from a 
group of particles. If we provisionally assume that the overall density distribution is characterized 
by a small number of discrete values that each represent a certain fraction of the total area, an 
assumption that comports with the qualitative sense of many investigators and that forms the 



foundation of the bimodal-distribution assumption of Colwell et al. (2006 20071 and Hedman et al. 
(2007), then it is reasonable for the variation we just described to transform the density histogram 
(Fig. [3J from a distribution containing only a small number of non-zero values to a distribution 
characterized by a small number of gaussians. 

Defining a single gaussian curve as y(x) = a§e~^ x ~ ax ^ 2 1 2 "^ , where x here is logr^n, the param- 
eters of each gaussian as stated in Fig. [3] are the height of the gaussian in units of the histogram 
(= ao), the location of the center of the gaussian in units of Td yn (= 10 ai ), and the full- width at 
half-maximum of logr^n (= 2\/2 In 2 • 02). The fit parameters are plotted in Fig. [4] as follows: 
the center location (= 10 ai+a 2/ 2 ) is now the characteristic value of Td yn for a log-normal distribu- 
tion, the gaussian height (= ao) is as in Fig. [3j and the integrated area under the gaussian curve 
(= \p2/K ■ aQd2) is normalized by the total area under the combined fitted curve. 

The close correspondence between the density histograms and our sum-of-three-gaussians mod- 
els justifies a posteriori our provisional assumption. Surprisingly, though, the density distribution 
is not bimodal but (except for the case with the lowest surface density) trimodal! Once the wakes 
are fully developed (i.e., for surface densities a > 30 g cm~ 2 in our simulations), the peak at high 
optical depth remains relatively stationary with a central value t ~ 2, and both the central value 
and the integrated area of the gaussian gradually increase with a (green in Fig. [4]). The largest 
fraction of the area is covered by very low optical depths, with fitted central values r ~ 0.04 that 
gradually decrease with increasing a (dark blue in Fig. [4]); these optical depths are so low that 
such regions are likely to be photometrically inactive in all cases. Optical depths closer to the T gap 
values inferred by previous authors characterize the middle peak in each histogram, with central 
values ranging from r ~ 0.1 to r ~ 0.25 (cyan in Fig. [4]). 



This result leads us to suggest re-interpreting the results of previous investigators (Colwell 
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Fig. 4. — Parameters of fitted gaussians from Fig. [3| Expressions are given in the text for (a) the 
center location, (b) the height, and (c) the integrated area normalized by the total area under the 
combined fit curve. The dotted line in panel (a) shows the input mean dynamical optical depth, 
calculated by dividing the sum of the cross-section areas of all particles by the total area. The 
colors of each curve correspond to the colors of the fitted gaussians in Fig. [3] Circles signify the 
simulations using a Borderies et al. (1984 1 coefficient of restitution law with v* = 0.001 cm s _1 , 



and triangles signify the simulation using a Bridges et al. (1984) law. 
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et al. 2006 2007 Hedman et al. 2007 Nicholson and Hedman 2009), whose models assume that a 



photometrically active T gap characterizes the entire area that is not occupied by the dense wakes. 



As we will show in Section 5.1 our results are likely consistent with this picture insofar as T gap is 



interpreted simply as the area- weighted average optical depth over the inter- wake regions; however, 
our results indicate that most of that area may be so sparsely populated that it is not photomet- 
rically active at all, while a smaller fraction of the inter-wake regions contains somewhat higher 
optical depths that cause the entire inter- wake region to average out to T gap - When one looks at 
movies of our self-gravity wake simulations, these regions of intermediate r most commonly arise 
from formerly dense wakes that are in the process of being disrupted. 

It is possible that the presence in our numerical simulations of large areas with very low optical 
depth might be a result of our use of a monodisperse size distribution (that is, we used a single 
particle size, namely R = 1.667 m). Preliminary work by Salo and Schmidt (2007) indicates that 
smaller particles may be more likely to escape the dense wake and spread into the inter- wake regions. 
Although we do not in this paper investigate simulations with a particle size distribution, this could 
be done easily enough by separating simulated particles into logarithmic size bins, perhaps with 
widths of a decade or a half-decade, assigning an average size to the particles within each bin, 
proceeding with a separate density estimation for each size bin, and finally combining the resulting 
estimated densities. On the other hand, there is considerable evidence that the clearing of small 
particles out of the inter- wake regions is in fact required in order to match brightnesses derived from 
TV-body wake simulations with the amplitude of the azimuthal brightness asymmetry observed in 



the A ring in both ground-based and spacecraft imaging data, as first suggested by Porco et al. 



(2003) and seen in the simulations of Porco et al. (2008). The clearing out of inter- wake regions 



occurs in models using a lossier ( "squishier" ) coefficient of restitution, which Porco et al.| (2008) 
argue are more realistic and provide a better fit to observations. Indeed, in our Figs. [| and g it 



is clear that the simulation using the "bounder" Bridges et al. (1984) coefficient of restitution law 



yields a much smaller peak (i.e., less area covered) at very low optical depths. 

As a check on the accuracy of our parameterization methods, in Fig. [5] we compare several 
methods of calculating the total optical depth of a simulated ring patch. The dotted line is the input 
mean dynamical optical depth, which would be equal to the actual mean photometric optical depth 
if particles were randomly distributed. The open symbols connected by a dashed line indicate the 
actual mean photometric optical depth, which can be easily calculated for the entire patch (though 
not locally); the transparency T is simply the fractional area not blocked by particles, measured 
by projecting particle positions onto a plane perpendicular to the line of sight, and Tphot = ~ 

InT. 



We note in passing that our calculated actual mean optical depth for the Bridges et al. simulation 



is 0.34; in Fig. 16 of Salo et al. (2004), whose simulation we intentionally attempted to replicate 
for purposes of comparison, those authors found a photometric optical depth of 0.36 for B = 90°. 

Ideally, the calculated total optical depths using our density estimation would be equal to the 
actual mean photometric optical depth. The solid symbols connected by a solid line indicate the 
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Fig. 5. — Solid symbols connected by solid line are parameterized mean optical depths, calculated 
by combining the central optical depths of each gaussian peak (Fig. [4^,) weighted by their relative 
integrated areas (Fig. |4j;). Open symbols connected by dashed line are actual mean photometric 
optical depth, calculated by dividing the area blocked by particles by the total area. Dotted line 
shows the input mean dynamical optical depth, calculated by dividing the sum of the cross-section 



areas of all particles by the total area. In all cases, circles signify the simulations using a Borderies 
et al. ( 1984 1 coefficient of restitution law with v* = 0.001 cm s _1 , and triangles signify the simulation 



using a Bridges et al. (1984) law. 



total optical depth obtained by combining^ the central locations of the gaussian peaks (Fig. |4p,) 
weighted by the integrated area under each gaussian curve (Fig. [4)3). If this final method of 
expressing the optical depth can be legitimized, then we will have succeeded in truly parameterizing 
our simulated self-gravity wakes, expressing their photometric properties by means of only a few 
numbers, whose dependence on surface density and other input factors can then be easily tracked. 

Firstly, we verify that the mean optical depth calculated from our parameterization is consistent 
with the mean optical depth calculated by simply adding together the continuum of density values 
weighted by the histograms shown in Fig. [3j These two are consistent within At ~ 0.003 (except 
for the first point, at a = 10 g cm -2 , for which At ~ 0.01), a small variation given that the curve 



Mil general, optical depths are properly combined by converting them to transparencies (T = e _T ), finding the mean of the 
transparencies weighted by their occurrence frequencies, and then converting the total transparency back into optical depth, 
thus: 

F 1 



' combined 



In 



N 



(i) 
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covers a range of At ~ 0.2. 

Secondly, however, we find that the parameterized mean optical depth (solid symbols) is slightly 
but uniformly too low compared to the actual mean optical depth (open symbols). The reason for 
this, we believe, is that we have heretofore used the dynamical optical depth Tdyn (which is easily 
calculated from dynamical simulations), rather than the photometric optical depth T p hot (which 
corresponds to observations). As we have mentioned, these two are roughly equal when particles 
are randomly distributed; this is because a gaussian probability of particle overlap (as seen when 
projected onto a plane perpendicular to any given line of sight) plays the same role in Td yn that the 
exponential plays in T p h ot . However, T p h ot > Td yn when the distance between particles is comparable 
to the particle size (i.e., at a high volume filling factor) because particles are now constrained as to 
the locations in space they can occupy. 



Salo and Karjalainen (2003) analyzed the photometric properties of simulated particle disks, 
finding that 

Tphot I 'Tdyn ^ 1 + kD, (2) 

where A: is a scalar of order unity and 

D = (4R/3H)T dyn (3) 

is the volume filling factor for vertical scale-height H and particle radius R. It is simple to obtain H 
from the results of our numerical simulations; it is twice the standard deviation of the z-coordinates 
of the particles, and we find it to range from 4 to 8 m. We attempted to calculate separate values 
of H for low-density and high-density bins but did not get a result that was statistically more 
meaningful than simply using a single value for each simulation; thus we have done the latter. All 
of our simulations used R = 1.667 m. 

The filling factors for each of our parameterized peaks are shown in Fig. [6^. Here we have 
used in place of Td yn the center location in optical depth of each peak (Fig. |4^). We then use the 
filling factor values to correct the peak locations via Eq. [2j as shown in Fig. [6)3. We show results 
for three values of k, namely 0, 0.5, and 1. The revised mean parameterized optical depth is shown 
in Fig. [6J:. Now we find that the mean optical depth given by our parameterization corresponds 
well with that calculated directly, especially for k = 1, similar to the value of k found to give the 



best results by Salo and Karjalainen (2003). 



The uncertainties in the parameter fits from which we obtained optical depth values (Figs. [3] 
and [4]) are difficult to quantify due to the high degree of correlation among the parameters. To 
estimate their order of magnitude, we fit two independent gaussians to the two most peak-like 
regions (namely, Td yn < 0.07 and Td yn > 1) of a characteristic histogram (namely, that with 
a = 45 g cm -2 ). We found, under the assumption that we had chosen the proper fitting function, 
that the parameters oq, «i, and a<i were robust within a few percent. 
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Fig. 6. — (a) Volume filling factor, D, calculated from Eq. [3] for each parameterized peak using 
values of R and H obtained for each simulation as a whole, and using the peak center location 
a s Tdyn- (b) Revised peak locations (compare Fig. [4^,) after applying Eq. |2j using k values of 0, 
0.5, and 1. (c) Revised mean parameterized optical depth (cf. Fig. [5|, calculated by combining 
the revised central optical depths of each gaussian peak (previous panel) weighted by their relative 
integrated areas (Fig. |4j;). 



5. Discussion 



We have now successfully parameterized our simulated self-gravity wakes, expressing their 
properties in terms of three optical depths that are weighted by the integrated area under each 
gaussian curve. After converting from dynamical to photometric optical depth by estimating the 
volume filling factor and using k = 1 for the scalar in Eq. [2j we find that the mean optical 
depth calculated through our parameterization corresponds well to that calculated directly from 
the simulations. We now proceed to apply this method to analysis of Cassini data from occultations 
and from imaging. 



5.1. Occultations 



In contrast to the bimodal optical-depth models employed by previous analysis of occultation 



data (Colwell et al. 2006, 2007 Hedman et al. 2007 Nicholson and Hedman 20091, we find for a 



ring patch in the mid-A ring with well- developed self-gravity wakes that slighly more than half of 
the area is taken up by space that for practical purposes is completely empty and photometrically 
inactive. The dense wakes themselves take up slightly less than one-fifth of the area, and the 
remaining one- fourth to two-fifths is characterized by the so-called T gap . Although simulations 
more detailed and accurate than those we performed might result in some quantitative changes 
to these results, the general idea of a trimodal distribution of optical depth in the rings deserves 
consideration. As we will show, it turns out that the two models perform nearly identically except 
at very low values of the ring opening angle B, which ranges from zero (seeing the rings edge-on) 
to 90° (seeing the rings face-on). 

A trimodal optical depth distribution should primarily affect the estimations of the gap optical 
depth Tg a p and the gap filling fractions F{ for the simplified wake models. To explore this issue, let 
us consider the transmission through the rings as a function of ring opening angle B in the case 
where the projection of the line of sight to the star onto the ringplane is parallel to the mean wake 



orientation, which was the geometry studied by Nicholson and Hedman (2009). In the formulation 



used by those authors, the fraction of the visible area taken up by inter-wake regions F 2 is equal to 
their parameter G/X, while Fi = 1 — F2 is ignored because T wa ^ e is taken to be infinite for practical 
purposes. Thus, for the simplified model of opaque wakes separated by gaps of constant optical 
depth (that is, the bimodal case) the transmission is given by 



T\B) = F 2 e 





I sinB 



(4) 



Similarly, for a model with a trimodal distribution of optical depths, the transmission is 

3 

T\B) = s £ j F i e- T *l^ B . (5) 



8=1 



where Fi and Tj are the area fractions and optical depths of the three components of the ring 
(e.g., Figs.]!]: andjGja, respectively). 
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Fig. 7. — Apparent normal optical depth r n as a function of the sine of the ring opening angle 
B for a bimodal (dashed) and a trimodal (solid) model. The bimodal model is calculated from 
Eq. [7] with F 2 = 0.8 and 

T gap — 0.09, while the trimodal model is calculated from Eqs. [5] and [6] 
with Fi = {0.2,0.3,0.5} and T gap = {4.0,0.2,0.02}, respectively. Specific values used for the area 
fractions Fi and optical depths Tj are shown in the plot. 




As discussed by Nicholson and Hedman (2009), it is useful to consider the apparent normal 
optical depth, 

T n = — sin B In T, (6) 
For the simplified (bimodal) wake model, this is a linear function of sinB: 

r h n = Tg ap - In F 2 sin B. (7) 

By contrast, the more complex trimodal model cannot be expressed so simply and will therefore 
yield a that is a nonlinear function of sin 13 . 

Apparent normal optical depths r n for both bimodal and trimodal models are plotted in Fig. [7J 
for values of the area fractions and optical depths roughly consistent with those found in Section |4| 
The bimodal model uses a gap fraction equal to the sum of the fractions of the two lower-optical- 
depth components in the trimodal model and a gap optical depth equal to the area-weighted mean 
optical depth in the two lower-optical-depth regions. For intermediate to moderately high ring 
opening angles, these two curves match reasonably well (up to a slight offset), which indicates that 
the optical depth of gaps derived from the occultation data is the mean optical depth of the regions 
outside the optically thick wakes. Some slight deviation from a straight line can be observed at 
the largest opening angles (sin£? ~ 1), as the occultations become marginally sensitive to the finite 
optical depth of the wakes. A more dramatic deviation from the linear trend, however, can also 
be seen at very low opening angles (sin I? < 0.2, which is to say B < 10°), where the moderate- 
optical-depth regions start to appear opaque, and thus the occultations begin to become capable 
of distinguishing between the very-low- and moderate-optical-depth components of the model. 
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Thus, occultations at very low opening angles may be able to discern variations in the optical 
depths within the gaps. Otherwise, existing analysis of occultation data should simply be inter- 
preted with Tg a p representing the area-weighted average optical depth within the gaps (or inter-wake 
regions), keeping in mind the possibility that there may be strong variations in optical depth within 
those inter- wake regions. 



5. 2. Imaging 



Porco et al. (2008) have already applied detailed direct photometric modeling, using a deter- 
ministic geometric ray-tracing method to shoot photons at simulated ring patches characterized by 
self-gravity wakes, to match the brightness of the ring as seen in Cassini images. Our semi- analytic 
method is not able to improve on their analysis in terms of quantitative properties, but we can 
speak to some trends that underlie the rings' photometric behavior. 

For images of the lit face of the rings, the large majority of the brightness comes from the dense 
wakes, which reflect more light than do regions of lesser optical depth. For images of the unlit face 
of the rings, on the other hand, most of the brightness comes from regions of intermediate optical 
depth, as the dense wakes are largely opaque and the nearly-empty regions do not scatter much 
light towards the camera. It would be instructive for future papers on direct photometric modeling 
to show at least a sample image of the brightness resulting from a ring patch; we predict that the 
brightness on the unlit side of the rings will be dominated by a relatively small area of intermediate 
optical depth, rather than being broadly spread over the inter-wake region. Furthermore, the trends 
discussed in this paper can help future direct photometric modeling to better tune their simulation 
parameters to match measured brightnesses. 

A more direct application of our parameterization is useful in investigating the peculiar photo- 
metric properties of so-called "propellers" . These local disturbances in the rings are thought to be 
due to the perturbing influence of embedded ~ 100-meter moonlets. Observed propellers (Tiscareno 



et al. 2006 2008 Sremcevic et al. 2007) generally show the disturbed region as brighter than the 



background ring brightness, though theoretical models ( Spahn and Sremcevic|200"0 Sremcevic et al. 



2002 Seifi et al. 2005 Lewis and Stewart 2009 ) indicate that the most prominent feature of the 



disturbance is a decrease in local density. For images of the unlit side of the rings, depending on the 
background optical depth, brightness sometimes decreases with increasing density as the rings be- 



come more opaque, but calculations of the photometric properties based on Chandrasekhar (1960) 



indicate that this is not the operating regime for the images in question, given their particular 



viewing geometry and the mean surface density as measured from spiral density waves (Tiscareno 



et al.|[2007 ). Furthermore, for images of the lit side of the rings, denser regions are always brighter 
than less dense regions. 



We apply our parameterization of self-gravity wakes to this problem by converting the three 
optical depth values for each simulation to brightnesses and then performing an area-weighted 
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Fig. 8. — Relative single-scattering brightness of a homogeneous slab of particles, calculated from 
Eqs. [8j (reflection, here shown as dashed line) and [9] (transmission, here shown as solid line) using 
(j, = 1 and fiQ = 0.4. 



average as before. Chandrasekhar ( |1960 ) derived for a homogeneous slab the single-scattering 



reflection (Rss) and transmission (T55): 

R SS oc 1 - e -^//WMo ( 8 ) 

and 

T ss oce- T/fl -e" T/w , (9) 

where [i and \i§ are the respective cosines of the emission angle (the angle of the direction towards 
the camera from the ring-plane normal, which is the complement of the ring opening angle B) and 
the solar incidence angle (the angle of the direction towards the Sun from the ring-plane normal) . 
We neglect the albedo and the phase function, which do not vary among the instances we consider 
and thus do not contribute to any trends. We will use fj, = 1, corresponding to an image in which 
the camera is looking straight down onto the face of the rings (which is to say B = 90°), because 
in carrying out our density-estimation method we projected ring particle positions onto the ring 
plane. To investigate behavior at other values of fi, one could redo the density estimation with ring 
particle positions projected onto a plane perpendicular to a different line of sight. Nevertheless, the 
general trends illuminated in our analysis here should be applicable to a broader range of viewing 
geometries. We use fio ~ 0.4, appropriate to the Sun's illumination of Saturn's rings in 2004 and 
2005, when the first images containing propellers were taken. The resulting values for Rss and Tss 
are plotted in Fig. [8} 



Our use of the single-scattering approximation is justified for images in reflected light (Dones 



et al. 1993) and for self-gravity wakes with their high filling factor (Salo and Karjalainen 2003) 



However, multiple scattering can significantly brighten images in transmitted light of non-wake 
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Fig. 9. — Calculated relative brightness for images in (a) reflection and (b) transmission, calculated 
from the optical depths given in Fig. [6)3 for k = 1 using the single-scattering brightnesses shown 
in Fig. [8] for a viewing geometry defined by n = 1 and ^0 = 0.4. Colors are for the three peaks, 
as in Figs. [3] and [4] while black is the combined brightness calculated by averaging the values for 
the three peaks weighted by the integrated areas given in Fig. |4j:. Dotted lines show the brightness 
associated with the input mean dynamical optical depth, which corresponds to the brightness a 
ring patch would have if self-gravity wakes were not present. 



regions if the ring reverts to a many-particle-thick structure (Salo and Karjalainen 20031. Thus, 
the dotted line in Fig. [9)3 may be even higher than we have shown, which only serves to strengthen 
the arguments made in the next section. 



5.3. Why propellers are bright 



The brightnesses resulting from our parameterized self-gravity wakes are shown in Fig. [9} We 
find for both reflection and transmission that the calculated brightnesses are quite low compared 
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to the brightnesses expected from an unperturbed ring patch lacking self-gravity wakes. In many 
cases, this is true even if the unperturbed patch has a lower input mean optical depth than the 
patch with self-gravity wakes. 

These results lend preliminary quantitative support to a mechanism previously proposed by 



Tiscareno et al. (2008), developing an idea first briefly suggested by Tiscareno et al. (2006), for why 



propellers appear as bright relative to the background ring. The argument in our previous paper 
was that self-gravity wakes tend to lock up ring material into a photometrically inactive state, and 
that propellers can release this material by locally and temporarily disrupting the wake structure. 



Thus, even though a propeller structure contains less material overall (Seifi et al. 2005; Lewis and 



Stewart 2009), it may contain more photometrically active material and thus be relative-bright. 



It should be noted that this hypothesis cannot be directly tested at this time by numerical 
simulations. To our knowledge, the most detailed simulations carried out to date on propellers 



are those by Lewis and Stewart (2009), but even their results do not fully incorporate the activity 



of self-gravity wakes, due to the great difficulty in simultaneously accounting for phenomena on 
very different lengthscales (a few meters for self-gravity wakes, several km for propellers). The 



numerical simulations and photometric models in the online supplement of Sremcevic et al. (2007) 



also are insufficient to either confirm or reject this hypothesis due to the severely limited extent of 
phase space explored. Investigating small-scale structure requires a large number of particles per 
unit area, while the large patch size required to investigate large-scale structure drives the total 
number of particles quite high. Our hope is that the semi-analytical treatment presented here will 
provide guidance to the highly computationally-intensive direct models that would be required to 
quantitatively address this question. 



Other possible hypotheses of course exist, such as the suggestion of Sremcevic et al. (2007) 



that densities within a propeller structure are enhanced by temporarily liberated ring-particle 
regolith, but it is not within the scope of this paper to study them. Our hypothesis was very briefly 
rebutted^ in the online supplement of Sremcevic et al. ( 2007[ ), who rejected the idea primarily 



because of a general sense that it should lead to "moonlet wakes"^ even brighter than the propeller 
gaps. However, many realistic simulations of propeller structures lack moonlet wakes (Lewis and 



Stewart 2009), and in any case, it may well be that both wakes and gaps are seen as bright in 



observed propellers. 

Our basic suggestion is that the continuum brightness of the ring is characterized by the solid 
black lines in Fig. |9j especially in the range of input mean surface densities from 40 to 45 g cm~ 2 



that typify the "propeller belt" region (Tiscareno et al. 2007), and that the brightness of well- 



developed propeller-shaped features is characterized by the dotted line. Although the dotted lines 



^in r esponse to informal di scussion during 2006, as well as the brief mention by 



See 



Tiscareno et al. 



Tiscareno et al. 



(2006) 



( 2008 ) and references therein for more in the context of propellers on the meaning of the term "moonlet 



wakes," originally due to |Showalter et al.| jl986| l, which has nothing to do with the "self-gravity wakes" discussed elsewhere in 



this paper, despite an unfortunate similarity in the names. 
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in Fig. [9] are above the solid black lines for any given surface density, the comparison must be 
made with the level of the dotted line at a lower surface density than for the solid black line, since 
the relative density inside well-developed propeller-shaped features is lower than the background - 
e.g., 10% to 30% in the models of Lewis and Stewart (20091. Propellers would be relative-bright 



features if the dotted lines at 10 to 15 g cm 2 have higher relative brightness than the solid black 
lines at the highest surface densities plotted. 



In keeping with the hypothesis of Tiscareno et al. (2006, 20081, Fig. 9p shows that all three 



parameterized optical depth values (colored solid lines) for images in transmission correspond to 
brightnesses lower than those of unperturbed rings with no wakes (dotted line) , even if the latter are 
considered to have significantly lower optical depths overall. The low values of the area-weighted 
average brightness make it quite likely that propellers may be seen as relative-bright features in 
images of the unlit side of the rings. For images in reflection (Fig. [9^,), the wakes themselves are 
quite bright (green symbols), and an average between only the dense wakes and the intermediate 
optical depth ("r gap ") would certainly result in wakes that are relative-dark features, but in our 
results the fact that a large fraction of the area is taken up by effectively empty space pushes the 
mean brightness (black solid line) downward to a regime in which propellers might be neutral or 
perhaps relative-bright features, when seen on the lit side of the rings. 

To be sure, we have not shown that this mechanism works in all cases. Factors that make 
it more likely to work include 1) higher densities within the propeller gaps, 2) higher background 
surface density a, and 3) lossier coefficient of restitution. In some scenarios, the propeller gaps might 
be brightened by this mechanism only enough to blend in with the background, or may remain 
relative-dark. It must also be remembered that our analysis is based on several approximations; 
thus, we cannot say in detail under which conditions this mechanism can or cannot work. Our 
purpose is to show that the mechanism is plausible, and to encourage more detailed work. 



These conclusions fit well with the observations. Tiscareno et al. (20081 found that propellers 



on both the lit and unlit faces of the rings were relative-bright features, and in fact found no 



relative-dark features in the "propeller belt" of the mid-A ring. Furthermore, Tiscareno et al. (2008) 



found no significant difference in propeller dimensions between lit-side (reflection) and unlit-side 
(transmission) images, indicating that the same basic structures were being seen in both geometries. 
Yet in morphology the observed propellers strongly resemble the regions of relatively low density 
seen in simulations, not the relatively dense "moonlet wakes" associated with the central moonlet. 



Preliminary observations of very large propellers by Tiscareno et al. (20091 do show, for the 



first time, some relative-dark regions within propellers. These "giant propellers" occur between the 
Encke Gap and the Keeler Gap, farther from the center of Saturn than the "propeller belt" in the 
mid-A ring that contains the smaller propellers observed previously. But even these observations 
can be incorporated into the above model simply as a matter of degrees, in effect becoming the 
exception that proves the rule. If relative densities inside the giant propellers are exceedingly 
low, then they will indeed be dark relative to the background in both reflection and transmission 
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(Fig. [9]) . Similarly, if relative densities inside the "moonlet wakes" are exceedingly high for the giant 
propellers, they may become so opaque that they are indeed relative-dark for images in transmission. 
However, neither of those circumstances affects our conclusion for the smaller propellers in the 



"propeller belt."** For these, we have found that propellers may very reasonably be found to be 
relative-bright features, both for images in reflection and for images in transmission, if they locally 
and temporarily disrupt the structure of self-gravity wakes in addition to decreasing the overall 
local density within the propeller-shaped structure. 



6. Summary and Conclusions 

We have developed a semi-analytic method of parameterizing simulations of self-gravity wakes 
in Saturn's rings, describing their photometric properties by means of only six numbers: three 
optical depths and three weighting factors. These numbers are obtained by fitting a sum of three 
gaussians to the results of a density-estimation procedure that finds the frequencies of various 
values of local density within a simulated ring patch. In order to account for the conversion from 
dynamical optical depth to photometric optical depth, we use the expression r p h ot / Td yn ~ 1 + kD 
( Salo and Karjalainen|[2003 ), where D is the volume filling factor (which we estimate); we find the 



best results by setting the scalar k = 1. 

Our first surprising result is that rings dominated by self-gravity wakes appear to be mostly 
empty space, with more than half of their area taken up by local optical depths around 0.01. Such 
regions will be photometrically inactive for all viewing geometries. While this result might be 
affected by our use of a monodisperse size distribution, we suggest that it might be robust due to 



the lower coefficient of restitution that we used in our iV-body simulations, following Porco et al 



(2008). 



If our models are in fact robust, then the bimodal density distribution assumed in the interpre- 
tation of occultation data by previous investigators should be replaced by a trimodal distribution. 
The practical results of this turn out to be minimal, as occultations can only distinguish between 
bimodal and trimodal models at very low opening angle. Existing analysis of occultation data 
should be interpreted with T gap representing the area-weighted average optical depth within the 
gaps (or inter- wake regions), keeping in mind the possibility that there may be strong variations in 
optical depth within those inter- wake regions. 

Applying our parameterization of self-gravity wakes lends preliminary quantitative support to 



the hypothesis of Tiscareno et al. (2008) that "propellers" observed in the mid-A ring are bright 
because of a local and temporary disruption of self-gravity wakes. Even though the overall local 
density is lower within the propeller-shaped structure surrounding an embedded central moonlet, 



"Also, background surface densities are lower beyond the Encke Gap than they are in the "propeller belt" (Spilker et al. 
2004 k 



- 25 - 



disruption of the wakes would flood these same regions with more photometrically active material, 
raising their apparent brightnesses in agreement with observations. We find that this mechanism 
can plausibly work for a wide range of input parameters. 

We hope that this hypothesis will eventually be directly tested by detailed numerical simula- 
tions. However, this is presently very difficult due to the high expense in terms of computational 
resources necessary to account for the small lengthscales appropriate for self-gravity wakes (a few 
meters) while simultaneously accounting for the very large lengthscales appropriate for propellers 
(several km). The former essentially requires a large number of particles per unit area, while 
the latter requires a large patch size. As computational resources increase to the point that such 
simulations become feasible, it will become possible to determine whether propellers really are 
characterized by a local and temporary disruption of self-gravity wakes, as well as whether the 
photometric properties of propellers can be explained using the hypothesis we have outlined. In 
the meantime, the results of our semi-analytic method give us reason to hope that the solution to 
the problem indeed lies in this direction. 
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